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ABSTRACT 

We  have  computed  theoretical  models  of  circumstellar  disks  for  the  classical  Be  stars  k  Dra,  /?  Psc,  and  v  Cyg. 
Models  were  constructed  using  a  non-LTE  radiative  transfer  code  developed  by  Sigut  &  Jones  (2007),  which  in¬ 
corporates  a  number  of  improvements  over  previous  treatments  of  the  disk  thermal  structure,  including  a  realistic 
chemical  composition.  Our  models  are  constrained  by  direct  comparison  with  long-baseline  optical  interferometric 
observations  of  the  Ha  -emitting  regions  and  by  contemporaneous  Ha  line  profiles.  Detailed  comparisons  of  our  pre¬ 
dictions  with  Ha  interferometry  and  spectroscopy  place  very  tight  constraints  on  the  density  distributions  for  these 
circumstellar  disks. 

Subject  headings:  circumstellar  matter  —  stars:  emission-line,  Be  —  stars:  individual  (k  Dra,  (3  Psc,  v  Cyg)  — 
techniques:  interferometric 
Online  material:  color  figures 


1.  INTRODUCTION 

Classical  Be  (or  B-emission)  stars  are  rapidly  rotating,  hot 
stars,  with  optical  spectra  that  show  hydrogen  emission  lines, 
and  frequently  emission  lines  from  singly  ionized  metals.  It  has 
been  recognized  since  the  days  of  Struve  (1931)  that  recombi¬ 
nation  in  a  flattened  disk  of  circumstellar  gas  can  reproduce  the 
basic  features  of  the  spectroscopically  observed  Ha  line  profiles 
of  Be  stars.  However,  it  is  now  generally  accepted  that  this  model 
is  too  simplistic.  The  detailed  mechanism(s)  that  creates  and  main¬ 
tains  Be  star  disks  remains  unclear  (Porter  &  Rivinius  2003; 
Owocki  2004),  although  rapid  rotation  of  the  central  B  star  cer¬ 
tainly  plays  a  role  (Townsend  et  al.  2004).  In  addition  to  this 
lack  of  a  successful  dynamical  model,  the  evolutionary  status 
of  Be  stars  and  related  stars  (such  as  B[e]  stars)  is  not  well  under¬ 
stood,  and  possible  evolutionary  connections  between  these  groups 
are  ambiguous.  These  problems  are  compounded  by  the  addi¬ 
tional  uncertainty  in  their  rotation  rates  (Cranmer  2005 ;  Townsend 
et  al.  2004;  Owocki  2004)  and  lack  of  understanding  of  how 
rotation  and  evolution  interplay.  The  suggestion  by  Townsend 
et  al.  (2004)  that  the  rotation  rates  of  Be  stars  may  be  systemati¬ 
cally  underestimated,  combined  with  new  interferometric  data, 
has  rekindled  interest  in  finding  suitable  dynamical  models  for 
these  stars. 

Until  recently,  the  bulk  of  the  observational  evidence  for  the 
presence  of  disks  surrounding  Be  stars  has  been  spectroscopic 
and  polarimetric  in  nature.  The  spectroscopic  evidence  includes 
the  profiles  of  the  emission  lines  observed  in  the  optical  and  1R 
spectra,  as  well  as  continuum  excesses  at  infrared  wavelengths. 
The  net  linear  polarization  observed  in  Be  stars  arises  via  elec¬ 
tron  scattering  in  the  nonspherical  distribution  of  gas,  as  first 
suggested  by  Coyne  &  Kruszewski  (1969).  However,  interfer¬ 
ometric  observations  that  spatially  resolve  the  disks  in  several 
wavelength  regimes  are  becoming  available  (see,  for  example, 
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Chesneau  et  al.  2005;  McAlister  et  al.  2005;  Tycner  et  al.  2006). 
Interferometric  observations,  in  concert  with  other  observables, 
allow  the  disk’s  physical  properties  to  be  determined  with  greater 
accuracy,  and  this  may  contribute  to  the  development  of  more 
successful  dynamical  models. 

Tycner  et  al.  (2006)  collected  long-baseline,  optical,  interfer¬ 
ometric  observations  of  the  classical  Be  stars  7  Cas  and  4>  Per  at 
the  Navy  Prototype  Optical  Interferometer  (NPOl)  and  inter¬ 
preted  these  observations  using  simple  models  of  the  intensity 
distribution  on  the  sky,  such  as  uniform  disks,  rings,  or  disks  with 
a  Gaussian  intensity  distribution.  In  this  paper,  we  extend  the 
work  of  Tycner  et  al.  (2006)  by  computing  detailed  theoretical 
models  for  the  intensity  distribution  on  the  sky  produced  by  the 
central  B  star  and  disk  of  a  Be  system.  We  then  compare  these 
models  to  observations  collected  at  NPOl. 

The  theoretical  models  are  computed  with  the  bedisk  code 
developed  by  Sigut  &  Jones  (2007).  Assuming  a  density  distri¬ 
bution  of  the  disk  gas,  the  code  can  compute  the  temperatures  in 
the  disk  given  the  energy  input  from  the  central  stars’  photo- 
ionizing  radiation  field  by  enforcing  radiative  equilibrium.  As 
the  disk  density  model  contains  several  adjustable  parameters, 
values  for  these  parameters  can  then  be  extracted  from  the  match 
to  observations.  We  present  new  observations  and  detailed  mod¬ 
els  for  the  Be  stars  k  Dra,  ti  Psc,  and  v  Cyg.  We  also  use  near- 
contemporaneous  observations  of  the  Ha  spectral  lines  for  these 
stars  as  additional  constraints  on  the  density  model.  As  we  shall 
demonstrate,  the  additional  spectroscopic  constraint  of  detailed 
line  profiles  is  often  an  important  ingredient  in  selecting  among 
models  consistent  with  interferometric  observations.  The  overall 
goal  of  this  study  is  to  place  tight  constraints  on  physical  condi¬ 
tions  in  the  circumstellar  regions  of  these  Be  stars. 

2.  THEORY 

The  theoretical  disk  models  presented  in  this  investigation 
were  constructed  using  the  bedisk  code  developed  by  Sigut  & 
Jones  (2007).  A  brief  overview  of  the  code  and  assumptions  rel¬ 
evant  to  this  work  are  presented  below.  The  reader  is  referred  to 
Sigut  &  Jones  (2007)  for  more  details. 

The  circumstellar  disk  is  assumed  to  be  axisymmetric  about 
the  star’s  rotational  axis  and  symmetric  about  the  midplane  of 
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the  disk.  In  this  case,  cylindrical  geometry  is  appropriate,  and 
we  use  R  for  the  radial  distance  from  the  star’s  rotation  axis  and 
Z  for  the  height  above  the  equatorial  plane.  The  radial  density 
distribution  in  the  equatorial  plane  is  assumed  to  be  given  by 
an  R~n  power  law,  following  the  works  of  Waters  (1986),  Cote 
&  Waters  (1987),  and  Waters  &  Cote  (1987).  At  each  radial 
distance  from  the  rotational  axis,  the  vertical  structure  of  the  disk 
is  determined  by  the  requirement  of  approximate  hydrostatic 
equilibrium  in  the  Z  direction,  which  balances  the  Z  component 
of  the  star’s  gravitational  acceleration  with  the  gradient  of  the  gas 
pressure  (see  Sigut  &  Jones  2007). 

Typical  values  for  the  radial  power-law  index  n  for  theoretical 
models  of  Be  stars  usually  fall  in  the  range  of  2-3.5,  based  on 
models  that  fit  the  1R  continuum  (Waters  1986).  Models  with  the 
highest  n  have  a  faster  decrease  in  the  density  distribution  with 
increasing  R.  The  density  distribution  within  the  disk  is  deter¬ 
mined  by  an  assumed  value  of  the  density  at  the  stellar  surface  in 
the  equatorial  plane,  po,  by  the  power-law  index  n,  and  by  the 
assumption  of  pressure  support  perpendicular  to  the  equatorial 
plane.  Since  po  and  n  significantly  affect  the  theoretical  predic¬ 
tions,  we  conduct  a  two-dimensional  parameter  search  by  varying 
these  input  parameters  over  all  reasonable  values  for  comparison 
with  interferometric  and  spectroscopic  observations. 

Given  this  density  model,  the  bedisk  code  is  used  to  find  the 
temperature  structure  of  the  disk  by  enforcing  radiative  equi¬ 
librium.  The  microscopic  heating  and  cooling  rates  implied  by 
a  gas  with  a  solar  chemical  composition  are  computed,  and  the 
temperature  that  balances  heating  and  cooling  is  found.  To  com¬ 
pute  the  atomic  level  populations  required  by  this  procedure,  the 
statistical  equilibrium  equations  are  solved  for  each  chemical 
element  included,  accounting  for  the  bound-bound  and  bound- 
free  collisional  and  radiative  processes  that  set  the  rates  in  and 
out  of  each  atomic  level. 

To  compute  the  hydrogen  line  profiles,  we  solved  the  transfer 
equation  along  a  series  of  rays  through  the  star-plus-disk  system 
as  viewed  at  an  angle  i  relative  to  the  line  of  sight  (where  i  =  0° 
is  pole-on,  and  i  =  90°  is  equator-on).  The  disk  was  assumed  to 
be  in  pure  Keplerian  rotation,  and  the  equatorial  rotational  ve¬ 
locity  of  the  star  was  chosen  so  that  the  measured  v  sin  i  of  the 
star  was  recovered  for  the  adopted  value  of  i.  For  rays  terminat¬ 
ing  on  the  stellar  surface,  we  adopted  a  photospheric  Ha  line 
profile  computed  using  the  synthe  code  and  the  Stark  broad¬ 
ening  routines  of  Barklem  &  Piskunov  (2003);  these  profiles  are 
based  on  the  hydrogen  populations  computed  for  the  appropriate 
LTE,  line-blanketed  model  atmosphere  adopted  from  Kurucz 
(1993).  For  rays  passing  through  the  disk,  the  equation  of  radiative 
transfer  was  solved  along  the  ray  using  the  short-characteristics 
method  of  Olson  &  Kunasz  (1987).  To  compute  the  Ha  opacity 
and  emissivity,  the  hydrogen  level  populations  computed  for  the 
thermal  solution  were  used.  Here,  hydrogen  was  represented  by  a 
15  level  atom,  and  all  implied  radiative  and  collisional  bound- 
bound  and  bound-free  processes  from  and  between  the  1 5  levels 
were  included.  The  radiative  bound-bound  rates  were  treated  by 
the  escape  probability  approximation,  in  which  the  net  radiative 
bracket  for  the  line  (Mihalas  1978)  is  replaced  by  a  single-flight 
escape  probability.  While  the  escape  probability  approximation 
can  provide  a  reasonable  description  of  the  line  thermalization, 
the  use  of  escape  probabilities  is  an  important  approximation.  The 
bedisk  code  assumes  static  escape  probabilities  based  on  the 
complete  redistribution  over  a  Doppler  profile  for  all  lines  (for  a 
justification  of  this  procedure,  see  Sigut  &  Jones  2007).  If  the 
functional  form  of  the  escape  probability  function  is  changed, 
for  example,  to  a  simple  r  dependence  on  optical  depth  for 
r»  1,  the  predicted  equivalent  width  of  Ha  can  change  by  ~2  A 


(due  to  the  change  in  the  thermal  structure  of  the  disk  and  in 
the  hydrogen  level  populations).  Thus,  in  matching  to  the  ob¬ 
served  Ha  profiles,  the  largest  uncertainty  is  likely  from  the 
uncertainty  in  the  theoretical  profile  and  not  from  the  obser¬ 
vational  uncertainties  (see  §  3  for  an  estimate  of  the  observational 
uncertainties). 

Our  approach  allows  the  full  dependence  of  the  Ha  line  emis¬ 
sivity  and  opacity  on  the  physical  conditions  throughout  the 
disk  to  be  retained.  The  Stark  broadening  routines  of  Barklem 
&  Piskunov  (2003)  were  also  used  to  compute  the  local  Ha  line 
profile  throughout  the  disk.  This  is  important,  as  these  routines 
can  handle  high  disk  densities  for  which  the  assumption  of  a 
simple  Gaussian  or  Voigt  profile  for  the  Ha  line  would  be  in¬ 
valid.  The  final  Ha  line  profile  from  the  unresolved  system 
was  obtained  by  summing  over  the  rays  weighted  by  their  pro¬ 
jected  area  on  the  sky.  The  Ha  line  profile  was  computed  over  a 
wavelength  bin  of  ±20  A  from  line  center  at  each  area  element 
on  the  sky.  To  obtain  the  final  Ha  line  profile  to  compare  with 
spectroscopic  observations,  the  profile  was  convolved  with  a 
Gaussian  of  FWHM  of  0.656  A  to  bring  the  resolving  power  of 
the  computed  profile  down  to  104  to  match  the  observations.  The 
predicted  interferometric  visibilities  also  followed  from  the  same 
numerical  model.  The  monochromatic  Ha  image  of  the  system 
projected  on  the  plane  of  the  sky  was  computed  by  integrating  the 
stellar  continuum  and  Ha  line  flux  over  a  1 50  A  spectral  window 
centered  at  Ha  to  represent  the  interferometric  spectral  channels 
(see  §  3).  Thus,  the  calculation  of  the  Ha  profiles  and  the  Ha  in¬ 
terferometric  visibilities  used  the  same  model  calculations. 

3.  OBSERVATIONS 

We  obtained  long-baseline  interferometric  and  spectroscopic 
observations  of  three  Be  stars:  k  Dra,  /3  Psc,  and  v  Cyg.  The 
interferometric  observations  were  acquired  using  NPOI  with 
baselines  ranging  from  18.9  to  64.4  m  in  length.  NPOI  is  de¬ 
scribed  in  detail  in  Armstrong  et  al.  (1998),  and  the  technique 
used  to  extract  interferometric  observables  from  the  spectral  chan¬ 
nel  containing  the  Ha  emission  line  has  been  discussed  by  Tycner 
et  al.  (2003).  Table  1  shows  all  the  interferometric  observations 
obtained  for  the  three  stars  and  the  number  of  data  points  ob¬ 
tained  on  each  night  in  the  spectral  channel  containing  the  Ha 
line.  The  interferometric  observations  of  all  stars  were  obtained 
over  a  period  of  less  than  a  month,  with  / 3  Psc  and  v  Cyg  observed 
during  2005,  and  k  Dra  observed  during  2006. 

To  complement  our  interferometric  observations,  we  have  also 
obtained  near-contemporaneous  spectroscopic  observations  of  all 
three  targets.  The  spectroscopic  observations  were  obtained  using 
an  Echelle  spectrograph  at  the  Lowell  Observatory’s  John.  S.  Hall 
telescope,  and  covered  the  spectral  region  around  the  Ha  line. 
The  properties  of  the  spectrograph  and  the  description  of  the  re¬ 
duction  methods  have  been  described  elsewhere  (see  §  3.2  in 
Tycner  et  al.  2006  and  references  therein).  The  continuum  signal- 
to-noise  ratio  near  Ha  is  more  than  200  in  the  spectra  of  all  three 
stars.  However,  the  equivalent  width  and  the  shape  of  the  line  is 
determined  with  respect  to  the  normalized  continuum  level.  Be¬ 
cause  the  line  profile  must  be  normalized  by  a  smoothly  varying 
function  that  is  fit  to  the  continuum  (or  where  one  thinks  the 
continuum  is  located),  the  choice  of  this  normalizing  function 
affects  the  shape  of  the  Ha  line.  In  addition,  the  presence  of 
telluric  lines,  as  well  as  their  temporal  variability  from  night  to 
night  (and  season  to  season)  can  affect  how  well  the  continuum 
level  is  determined;  we  expect  this  to  be  the  dominant  source 
of  uncertainty  associated  with  the  shape  of  the  observed  Ha 
line.  We  have  estimated  the  magnitude  of  this  effect  by  analyzing 
a  set  of  reduced  spectra  obtained  from  the  same  raw  spectrum, 
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TABLE  1 

Observing  Log  for  Interferometric  Observations 


UT  Date 

k,  Dra 

Number  of  Points 

/ 3  Psc  v  Cyg 

2005  Aug  28 . 

12 

8 

2005  Aug  29 . 

16 

12 

2005  Aug  30 . 

20 

20 

2005  Aug  31 . 

20 

20 

2005  Sep  1 . 

20 

20 

2005  Sep  5 . 

8 

4 

2005  Sep  13 . 

16 

16 

2005  Sep  14 . 

23 

20 

2005  Sep  15 . 

20 

20 

2005  Sep  16 . 

21 

24 

2005  Sep  17 . 

8 

9 

2005  Sep  26 . 

16 

28 

2006  Feb  25 . 

30 

2006  Feb  26 . 

6 

2006  Mar  5 . 

52 

2006  Mar  6 . 

24 

2006  Mar  14 . 

36 

2006  Mar  16 . 

76 

2006  Mar  19 . 

52 

Total: 

276 

200 

201 

but  all  with  slightly  different  quadratic  functions  fitted  to  slightly 
different  continuum  regions.  The  variations  in  the  resulting  nor¬ 
malized  continuum  were  at  the  2%— 3%  level;  to  be  conservative, 
we  adopt  a  3%  uncertainty  for  our  continuum  level  determina¬ 
tion  and,  in  turn,  the  shape  and  equivalent  width  measure  of  the 
observed  Ha  emission  line. 

For  /?  Psc  and  v  Cyg,  the  spectra  were  obtained  on  2005 
September  1 6,  and  for  n  Dra  the  spectrum  was  obtained  on  2006 
March  17;  therefore,  in  all  three  cases  the  spectroscopic  obser¬ 
vations  overlapped  the  interferometric  observing  runs  (Table  1). 
The  equivalent  widths  of  the  observed  Ha  spectra  lines  are 
—  15.3  A,  —27.3  A,  and  —22.2  A  for  (3  Psc,  v  Cyg,  and  k  Dra, 
respectively.  (We  adopt  the  usual  convention  that  a  negative  in 
the  quoted  value  of  the  equivalent  width  indicates  a  net  line 
flux  above  the  continuum.) 

4.  RESULTS 

There  is  considerable  uncertainty  in  the  assigned  spectral  types 
for  Be  stars  due  to  rapid  rotation  and  spectral  variability.  The 
effects  of  gravity  darkening  and  the  potential  obscuration  of  por¬ 
tions  of  the  stellar  surface  by  the  disk  further  compound  the  un¬ 
certainty  in  classification  and  the  corresponding  stellar  parameters. 
Since  the  stellar  parameters  fix  the  photoionizing  radiation  field 
that  is  assumed  to  be  the  sole  source  of  energy  input  into  the 
circumstellar  disk,  realistic  stellar  parameters  are  essential  to 
the  construction  of  reasonable  disk  models.  For  each  of  the  stars 
in  this  study,  we  searched  the  literature  in  order  to  find  the  most 


reasonable  stellar  parameters  available.  See  Table  2  for  the  adopted 
stellar  parameters  for  k  Dra,  6  Psc,  and  v  Cyg.  The  rationale  for 
adopting  these  particular  stellar  parameters  is  explained  and  com¬ 
pared  with  other  published  values  below. 

4.1.  n  Draconis 

k  Dra  is  a  known  spectroscopic  binary  and  giant  star  of  in¬ 
termediate  spectral  type.  It  is  variable  on  a  variety  of  timescales 
(see  Hirata  1995  for  a  detailed  discussion  of  these  variations). 
Saad  et  al.  (2004)  modeled  the  central  star  of  k  Dra  by  com¬ 
paring  a  variety  of  spectroscopic  and  photometric  observations 
to  a  grid  of  non-LTE  models.  We  have  adopted  their  derived 
best-fit  stellar  parameters  for  our  study.  We  used  their  param¬ 
eters  of  TeS,  log  g,  and  mass  to  calculate  the  stellar  radius.  Gies 
et  al.  (2007)  also  used  these  parameters  for  the  interpretation  of 
ik'-band  observations  of  k  Dra  collected  at  the  CHARA  Array 
interferometer.  The  parameters  of  Saad  et  al.  (2004)  are  also  in 
reasonable  agreement  with  stellar  parameters  listed  for  the  spec¬ 
tral  class  B5  111  in  Cox  (2000). 

Figure  1  shows  the  density  parameter  grid  for  the  disk  of 
k  Dra.  This  grid  was  selected  after  running  a  series  of  models 
over  a  much  wider  range  in  po  and  n  on  a  coarser  grid.  A  com¬ 
parison  of  the  observed  and  predicted  Ha  equivalent  widths 
allowed  us  to  determine  a  finer  grid  to  sample  for  detailed  com¬ 
parisons  with  interferometric  observations.  The  final  density  pa¬ 
rameter  grid  consisted  of  278  models.  We  have  compared  the 
output  from  each  model  to  the  interferometric  observations  using 
the  technique  described  in  detail  by  Tycner  et  al.  (2008).  The 
reduced  x2  values  corresponding  to  these  models  are  presented 
in  Figure  2  as  a  function  of  po  (g  cnr3)  and  n.  The  darkest 
patches  in  this  figure  represent  the  models  with  reduced  X2  <  2. 
The  white  plus  signs  correspond  to  models  with  the  lowest  val¬ 
ues  (ranging  from  1.18  to  1.29)  of  the  interferometric  reduced 
X2-  In  total,  there  are  25  models  with  reduced  x2  <  2  for  this  star; 
therefore,  it  is  necessary  to  further  constrain  our  models  in  order 
to  determine  the  best-fit  model. 

Figure  3  shows  the  predicted  Ha  equivalent  widths  as  a  func¬ 
tion  of  log  po,  and  the  power-law  index  n  for  the  entire  grid  of 
models.  The  predicted  Ha  profiles  were  constructed  assuming 
an  inclination  of  i  —  35°.  We  note  that  this  value  of  i  is  consis¬ 
tent  with  other  values  determined  from  spectroscopy  and  polar¬ 
ization  analyzes  presented  in  the  literature.  For  example,  Clarke 
(1990)  estimates  i  =  23°  based  on  the  observed  polarization  an¬ 
gle,  and  Juza  et  al.  (1991)  find  i  =  35°-45°  based  on  spectros¬ 
copy.  We  also  experimented  by  varying  i  and  found  that  an 
inclination  angle  for  this  system  of  35°  seemed  to  produce  pro¬ 
files  that  matched  the  observed  Ha  line.  We  also  note  that  al¬ 
though  changing  /  by  ±  1 0°  changes  the  shape  of  the  Ha  profile, 
the  effect  on  the  equivalent  width  is  <  1  A.  The  lightest  colored 
symbols  are  closest  to  the  value  of  the  observed  equivalent  width 
of  —22.2  A.  There  are  14  models  that  predict  an  Ha  equivalent 
width  within  ±2  A  of  the  observation,  and  these  are  repre¬ 
sented  by  white  circles.  Models  that  correspond  to  predicted  Ha 


TABLE  2 

Adopted  Stellar  Parameters 


Name  of  Star 

HR  Number 

log9 

Spectral  Type 

Stellar  Radius 
(Re) 

Stellar  Mass 
(M0) 

Teff 

(K) 

v  sin  i 
(km  s"1) 

Reference 

k,  Dra . 

HR  4787 

3.5 

B6  Illpe 

6.4 

4.8 

14000 

170 

i 

/3  Psc . 

HR  8773 

4.0 

B6  Ve 

3.6 

4.7 

15500 

90  ±  15 

2 

v  Cyg . 

HR  8146 

4.0 

B2  Vne 

4.7 

6.8 

19800 

173  ±  10 

3 

References. — (1)  Saad  et  al.  2004;  (2)  Levenhagen  &  Leister  2004;  (3)  Neiner  et  al.  2005. 
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Fig.  1 . — Density  parameter  grid  for  k  Dra  as  a  function  of  po  (g  cm  3)  and  n. 
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Fig.  3. — Predicted  Hck  equivalent  widths  as  a  function  of  log  po  and  the  power- 
law  index,  n,  for  the  entire  grid  of  models.  The  lightest  colored  symbols  are  closest 
to  the  value  of  the  observed  equivalent  width  of —22.2  A.  The  models  that  predict 
an  Ha  equivalent  width  within  ±2  A  of  the  observation  are  represented  by  white 
circles.  The  models  that  correspond  to  predicted  Ha  emission  greater  than  the  ob¬ 
servation  are  indicated  by  plus  signs  within  the  symbol. 


equivalent  widths  greater  than  the  observed  emission  have  a  plus 
sign  within  the  corresponding  symbol  in  Figure  3. 

We  note  that  this  subset  of  14  models  does  not  represent  the 
14  best  models  using  interferometry  alone.  Although  this  set 
matches  the  observed  Ha  equivalent  width  within  ±2  A,  if  we 
use  interferometric  observations  as  a  constraint,  this  same  set 
corresponds  to  a  reduced  x2  from  1.18  (the  minimum  value)  up 
to  16.0.  In  fact,  6  of  the  14  models  in  this  set  had  a  reduced  x2 


2  3  4  5 


n 

Fig.  2. — The  \2  values  for  computed  models  as  a  function  ofn  and  po  (gem  '3) 
for  k  Dra.  The  contour  plot  shows  the  reduced  \2  values  in  steps  of  1 .0.  Only  a 
region  of  reduced  \2  <  6  is  shown.  The  darkest  regions  on  the  figure  represent 
models  with  reduced  \2  <2.  The  four  white  plus  signs  correspond  to  models 
with  the  lowest  values  of  reduced  \2 ,  as  determined  from  interferometry.  The 
diamonds  correspond  to  models  that  have  predicted  Ho  equivalent  widths  ±2  A 
of  the  observation. 


greater  than  2.0.  It  is  clear  from  a  comparison  of  Figure  2  and 
Figure  3  that  the  results  from  Hct  interferometry  and  Her  line  pro¬ 
file  are  complementary,  but  that  each  provides  extra  information 
to  help  reduce  the  degeneracy  in  finding  the  best-fit  model. 

Multiple  pairs  of  n  and  po  that  produce  Ha  profiles  with  equiv¬ 
alent  widths  that  match  the  observations  within  uncertainties  can 
be  found.  However,  not  all  of  these  are  the  best  fits  in  terms  of 
interferometry.  That  is,  it  is  possible  to  match  the  equivalent 
width  for  a  range  of  n  by  an  appropriate  adjustment  in  po,  but 
these  combinations  do  not  necessarily  have  the  correct  density 
distribution  as  a  function  of  radial  distance  from  the  star.  How¬ 
ever,  by  combining  results  from  spectroscopy  and  interferometry, 
we  can  find  the  best  pairs  of  n  and  po  for  k  Dra. 

Next,  we  compared  the  subset  of  eight  models  that  represent 
the  best  fits  from  both  interferometry  and  Ha  equivalent  width. 
A  subset  (four  models  with  the  lowest  reduced  y2  from  interfer¬ 
ometry)  of  the  corresponding  Ha  profiles  is  shown  in  Figure  4. 
Note  the  variation  in  the  shape  of  the  profile  as  n  is  increased 
from  2.7  in  the  leftmost  panel  to  4.2  in  the  rightmost  panel.  The 
predicted  equivalent  widths  are  —21.9,  —23.5,  —22.4,  and  —20.3  A, 
with  corresponding  y2  values  from  interferometry  of  1.18,  1 .29, 
1.18,  and  1.18,  respectively.  The  model  with  n  =  2.1  does  not  have 
the  correct  line  profile  shape  compared  to  the  observed  profile. 
There  is  too  much  absorption  at  line  center,  and  the  profile  is 
clearly  too  narrow  in  the  wings.  All  of  the  predicted  profiles  in 
Figure  4  have  an  absorption  feature  at  line  center.  This  feature  is 
most  pronounced  for  low  n  and  gradually  decreases  for  models 
with  higher  n.  If  we  increase  the  inclination,  this  absorption  fea¬ 
ture  becomes  more  and  more  pronounced  as  the  system  is  viewed 
closer  and  closer  to  edge-on.  If  we  decrease  i,  so  that  the  system 
is  viewed  close  to  pole-on,  the  absorption  feature  remains  for  the 
models  with  the  lowest  n.  In  addition,  with  n  =  2.7,  the  wings  of 
the  predicted  profile  are  not  broad  enough  for  any  reasonable  val¬ 
ues  of  i.  The  panel  with  the  model  corresponding  to  n  =  3.8  and 
n  =  4.0  in  Figure  4  shows  that  the  predicted  line  is  too  strong, 
which  results  in  an  Ha  equivalent  width  that  is  too  large.  Mod¬ 
els  corresponding  to  «  >  4.2  have  profiles  with  wings  wider  than 
the  observation.  The  model  that  corresponds  to  the  profile  with 
n  =  4.2  matches  the  shape  of  the  observed  profile  best,  and  has  a 
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Fig.  4. — Observed  Hen  line  for  n  Dra  compared  with  4  models  with  different  values  of  n  and  po.  These  models  correspond  to  the  best  fit  from  both  interferometry  and 
Ha  modeling.  Solid  lines  represent  the  model  profiles,  and  the  circles  represent  the  observed  line.  All  model  profiles  were  obtained  with  a  disk  inclination  of  35°.  The 
reduced  \2  values  from  interferometry  for  these  models  are  1.18, 1 .29,  1.18,  and  1.18,  from  left  to  right.  [See  the  electronic  edition  of  the  Journal for  a  color  version  of 
this  figure .] 


minimum  reduced  x2  of  1 . 1 8  from  interferometry.  We  prefer  this 
model  with  n  =  4.2  and  po  =  1.5  x  10~1()  g  cnT3  in  terms  of 
both  spectroscopy  and  interferometry.  Figure  5  shows  the  inter¬ 
ferometric  observations  compared  with  this  model  for  t;  Dra. 

We  note  that  as  previously  discussed,  the  uncertainty  due  to 
the  approximations  in  the  computational  code  (see  §  2)  is  larger 
than  that  in  the  observed  Ha  profiles  (see  §  3).  We  searched  the 
literature  to  find  the  best  possible  stellar  parameters  for  the  disk 
systems  studied  in  this  investigation  (see  §§  4. 1-4.3).  However, 
the  stellar  parameters  for  Be  stars  are  generally  not  well  estab¬ 
lished,  and  for  some  stars  even  the  spectral  type  may  not  be  ac¬ 
curate.  We  now  wish  to  test  the  dependence  of  model  parameters 
on  our  results.  We  experimented  by  varying  the  stellar  parame¬ 
ters  in  turn,  to  see  how  these  changes  affect  the  predicted  Ha; 
equivalent  widths  and  the  interferometric  fits  for  the  star  k  Dra. 
Similar  types  of  uncertainties  would  be  expected  for  the  other 
two  systems  presented  in  this  paper.  One  might  expect  that  the 
Teff  of  the  central  star,  which  supplies  the  disk  with  energy, 
would  be  the  source  of  the  greatest  uncertainty.  We  changed  the 
TeS  by  ±1000  K  and  reproduced  all  of  the  profiles.  This  resulted 
in  a  change  in  the  equivalent  width  of  the  Ha  line  by  a  maximum 
of  ~2  A.  To  assess  the  effect  of  uncertainties  when  the  funda¬ 
mental  stellar  parameters  were  varied  in  combination,  we  per¬ 
formed  a  Monte  Carlo  simulation  in  which  new  stellar  parameters 
were  randomly  realized.  These  parameters  were  used  to  compute 
a  new  circumstellar  disk  model,  corresponding  to  our  preferred 
model  of  n  =  4.2  and  po  =  1.5  x  10~10  g  cm~3.  This  procedure 
was  repeated  to  generate  new  Ha  profiles  and  equivalent  widths. 


Spatial  Frequency  (106  cycles/radian) 

Fig.  5. — Interferometric  data  and  model  with  n  =  4.2  and  po  =  1.5x10  10  g 
cm-3  for  k  Dra.  The  dashed  line  represents  the  central  star,  and  the  solid  line  is 
obtained  by  taking  a  Fourier  transform  of  a  synthetic  image  of  the  preferred  model. 
[See  the  electronic  edition  of  the  Journal  for  a  color  version  of  this  figure.] 
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Errors  assumed  for  the  stellar  parameters  were  ±25%  in  mass, 
±25%  in  radius,  and  ±  1 0.7%  in  reff  (corresponding  to  ±  1 500  K); 
each  stellar  parameter  was  assumed  to  be  equally  probable  within 
its  errors.  Using  a  Gaussian  distribution  for  the  errors,  we  found 
that  our  results  did  not  change  significantly.  The  mean  of  the  Ha 
equivalent  width  was  —20  A,  with  a  standard  deviation  of  2.65  A. 

We  have  one  star  in  common  with  the  work  of  Gies  et  al. 
(2007):  k  Dra.  They  use  a  disk  density  model  very  similar  to 
ours  and  derive  two  sets  of  parameters  («,  po )  based  on  fits  to 
W'-band  interferometric  visibilities:  n  =  0.16  ±  0.43,  po  =  2.8  x 
10-13  g  cnr3  for  a  single-star  model,  and  n  =  0.67  ±  0.36, 
po  =  6.2  x  10-13  g  cnr3  for  a  binary  model.  Gies  et  al.  (2007) 
prefer  the  binary  fit  to  the  single-star  model  because  of  the  pe¬ 
culiar  solution  they  obtain  in  the  single-star  case;  namely,  a 
very  small  exponent  in  the  radial  density  power  law  (n  =  0.16) 
and  a  best  solution  with  i  =  72°  ±  18°,  which  is  at  odds  with 
the  resolved  Ha  line  profile.  Nevertheless,  the  particular  prop¬ 
erties  of  the  introduced  binary  are  not  constrained  by  any  in¬ 
dependent  observations.  Both  of  these  solutions  seem  at  odds 
with  our  results.  It  is  well  documented  that  the  Ha  equivalent 
width  and  the  infrared  excess  are  correlated  (see,  for  example, 
van  Kerkwijk  et  al.  1995).  Saad  et  al.  (2004)  present  a  figure 
that  shows  the  change  in  Ha  equivalent  width  from  the  early 
1970s  to  2004  for  n  Dra  (see  Saad  et  al.  2004,  Fig.  Ad).  Al¬ 
though  there  is  certainly  a  scatter  of  data  values  in  their  figure, 
one  can  easily  see  from  this  figure  that  the  Ha  equivalent  width 
varies  by  ~  1 5  A  ( from  about  —5  to  about  —20  A)  in  a  period  of 
~6  yr.  This  is  a  substantial  change  in  the  Ha  emission  line,  and 
we  expect  that  there  would  be  a  corresponding  increase  in  the 
infrared  excess.  The  interferometric  observations  for  k  Dra  pre¬ 
sented  in  Gies  et  al.  (2007)  were  obtained  in  2005  April  (3  dates) 
and  December  (2  dates).  Their  spectroscopic  observations  were 
obtained  in  2005  April.  Our  contemporaneous  observations  were 
obtained  in  2006.  The  difference  between  our  results  and  those 
of  Gies  et  al.  (2007)  is  partially  accounted  for  by  the  fact  that  we 
model  this  variable  star  at  different  times.  This  further  supports 
the  requirement  of  contemporaneous  observations.  Other  pos¬ 
sible  causes  for  the  difference  in  results  could  be  due  to  the  fact 
that  the  Ha  observations  are  sensitive  to  different  regions  within 
the  disk,  as  well  as  the  different  spatial  resolutions  of  the  two 
interferometers. 

There  are  also  other  differences  between  our  calculations  and 
those  of  Gies  et  al.  (2007)  that  should  be  highlighted.  As  noted 
in  §  2,  the  calculation  of  the  Ha  visibilities  and  the  Ha  line 
profiles  in  our  work  are  obtained  from  essentially  the  same  cal¬ 
culation  that  determined  the  radiative  equilibrium  temperature 
structure  for  the  disk,  given  the  adopted  density  model.  The  Ha 
opacity  and  emissivity  at  each  point  in  the  disk  were  computed 
from  the  local  temperature,  pressure,  and  radiative  field  by  solving 
the  set  of  statistical  equilibrium  equations  for  a  15  level  hydro¬ 
gen  atom.  While  our  fits  to  the  observed  Ha  line  profiles  are  not 
perfect,  they  are  predicted  naturally  from  the  models  and  the  as¬ 
sumption  of  pure  Keplerian  rotation  by  the  disk.  To  fit  the  Ha , 
Hy,  and  Bra  lines,  Gies  et  al.  (2007)  found  it  necessary  to  con¬ 
volve  their  profiles  with  a  unit-normalized  Lorentzian  to  broaden 
the  line  wings;  they  discuss  numerous  possible  sources  for  this 
extra  broadening.  In  the  current  work,  we  convolved  our  line 
profiles  only  with  a  unit-normalized  Gaussian  with  a  FWHM  of 
0.656  A  to  bring  the  resolving  power  of  the  computed  profile  to 
104  to  match  the  observations.  Among  the  reasons  cited  by  Gies 
et  al.  (2007)  as  possible  explanations  for  broad  Ha  wings,  they 
discuss  the  possibility  of  Stark  broadening  for  the  hydrogen  lines, 
but  then  note  that  the  densities  they  find  are  too  low  for  Stark 
broadening  to  be  important.  This  is  consistent  with  their  best- 


fit  densities,  which  are  about  1-2  orders  of  magnitude  smaller 
than  those  derived  here. 

Previous  studies  have  shown  that  self-consistent  disk  thermal 
structures  are  crucial  in  order  to  interpret  observations  correctly 
(see,  for  example,  Carciofi  &  Bjorkman  2006;  Sigut  &  Jones 
2007).  We  computed  the  thermal  structure  for  k  Dra  with  our 
bedisk  code  using  the  model  parameters  n  =  0.67  ±  0.36,  po  = 
6.2  x  10~13  g  cm~3,  andTG  =  67  R@  listed  in  Gies  et  al.  (2007). 
We  find  substantial  variation  in  disk  temperature  near  the  star 
and  equatorial  plane,  and  an  overall  average  disk  temperature  of 
~  13, 000  K.  Gies  et  al.  (2007)  adopt  an  isothermal  disk  temper¬ 
ature  of  ~9300  K  for  k  Dra  based  on  previous  studies  (Carciofi 
&  Bjorkman  2006),  but  we  note  that  neither  this  study  nor  Sigut 
&  Jones  (2007)  investigates  disk  densities  as  low  as  the  Gies  et  al. 
(2007)  best-fit  model.  Finally,  we  note  that  our  value  of  n  =  4.2 
for  k  Dra  is  in  agreement  with  typical  values  previously  predicted 
in  the  literature  for  Be  star  disks  (see,  for  example.  Waters  1986; 
Tycner  et  al.  2008). 

4.2.  (3  Pisces 

f3  Psc  is  a  B6  Ve  star  that  has  long  been  known  to  exhibit 
emission  in  the  hydrogen  lines  (see,  for  example,  Merrill  et  al. 
1925).  This  star  is  similar  in  spectral  type  to  k  Dra,  but  it  is  a 
main-sequence  dwarf.  We  have  included  f3  Psc  in  this  investi¬ 
gation  because  comparing  the  disk  models  with  interferometric 
observations  led  to  unique  density  parameters,  and  it  is  useful  to 
compare  the  results  for  this  star  with  those  for  k  Dra.  We  adopt 
the  stellar  parameters  from  Levenhagen  &  Leister  (2004)  and 
note  that  these  parameters  are  also  consistent  with  interpolated 
values  in  Cox  (2000).  The  interferometric  observations  are  from 
12  nights  between  2005  August  28  and  September  26  (see 
Table  1);  the  Ha  spectroscopy  for  (3  Psc  (and  v  Cyg)  were 
obtained  on  2005  September  16  and  therefore  are  definitely  rep¬ 
resentative  of  the  Ha  line  during  the  period  of  the  interferomet¬ 
ric  observations.  We  note  that  a  subsequent  observation  of  the 
Ha  line  for  (3  Psc  on  2006  October  3  had  a  20%  weaker  equiv¬ 
alent  width. 

The  procedure  for  selecting  the  grid  of  models  to  be  compared 
in  detail  with  interferometric  observations  was  obtained  by  the 
same  selection  process  used  for  k  Dra.  The  final  parameters  to  be 
searched  covered  the  range  n  =  1.5-5. 5  and  po  =  1.0  x  10~12- 
8.0  x  10~9  g  cm~3.  In  this  range,  193  models  were  constructed 
for  comparison  with  interferometric  observations.  The  correspond¬ 
ing  x2  values  had  a  minimum  of  0.83,  corresponding  to  the 
model  with  n  =  2.2  and  po  =  4.0  x  10~12  g  cm~3,  and  a  max¬ 
imum  of  438,  corresponding  the  model  with  n  =  1 .7  and  po  = 
1.0  xl0~n  g  cm~3.  The  x2  values  revealed  only  one  other 
model,  n  =  2.1  andpo  =  3.0  x  10  I2,  which  was  statistically  the 
same  as  the  best  fit.  For  this  model,  x2  =  0.84. 

In  order  to  check  the  predicted  Ha  profiles,  we  had  to  adopt  a 
value  for  i.  The  observed  sharp  emission  line  profile  (shown  in 
Fig.  6)  suggests  a  small  i.  However,  this  star  has  a  v  sin  i  of  90° 
km  s_1  (see  Table  2),  so  for  very  small  i,  the  star  will  exceed  its 
critical  velocity.  The  estimated  critical  velocity  of  a  star  of  this 
spectral  type  is  ~416  km  s_1  (Porter  1996).  Therefore,  we  argue 
that  the  inclination  of  this  system  must  be  at  least  15°,  and  we 
have  adopted  a  value  of  i  =  20° . 

We  calculated  the  Ha  equivalent  widths  for  a  subset  (~20 
models)  chosen  by  the  best-fit  x2  values,  and  we  considered  x2 
values  up  to  ^  1 .52.  Statistically,  the  two  models  corresponding 
to  x2  values  of  0.83  and  0.84  are  equally  good,  while  the  other 
models  in  the  subset  with  larger  values  of  x2  are  poorer  fits  to  the 
interferometric  data.  Nevertheless,  to  be  complete,  we  computed 
profiles  for  every  model  in  the  subset.  The  equivalent  widths  for 
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Fig.  6. — Observed  Ha  line  for  (3  Psc  compared  with  model  predictions  for  a  series  of  models  showing  the  changes  in  the  profile  as  n  is  increased.  Solid  lines  represent 
the  model  profiles,  and  the  circles  represent  the  observation.  [See  the  electronic  edition  of  the  Journal  for  a  color  version  of  this  figure .] 


the  two  models  corresponding  to  the  minimum  \2  values  had 
reasonable  line  shapes  compared  to  the  observed  profile,  but 
the  model  corresponding  to  n  =  2.2  had  an  equivalent  width  of 
—25.0  A,  while  the  model  corresponding  to  n  =  2.1  had  an 
equivalent  width  of  —2 1 .6  A,  which  is  closer  to  the  value  for  the 
observed  spectra  of  —15.3  A. 

Figure  6  {left)  shows  the  Ha  profile  for  the  best-fit  model 
(n  =  2.1,  po  =  3.0  x  10-12  g  cm~3)  in  terms  of  interferometric 
observations  and  spectroscopy.  Models  corresponding  to  other 
values  of  n  did  not  produce  Ha  profiles  that  agreed  with  obser¬ 
vations.  Figure  6  shows  the  behavior  of  the  Ha  profile  with  in¬ 
creasing  The  reader  can  clearly  see  that  the  wings  of  the  line 
in  the  middle  and  right  panels  are  too  wide,  resulting  in  a  pre¬ 
dicted  emission  that  is  too  large.  These  profiles  correspond  to 
models  with(«  =  3.6,  po  =  1.5  x  10-10  gem-3,  equivalent  width 
of  —30.4  A)  and  (n  =  4.0,  po  =  4.0  x  10-10  g  cm-3,  equivalent 
width  of  —23.4  A),  and  have  interferometric  reduced  x2  values 
of  1 . 1 3  and  1.21,  respectively.  Models  with  n  <2.1  had  large  %2 
values  based  on  the  interferometric  results,  and  were  eliminated. 
Figure  7  shows  the  interferometric  observations  for  our  best-fit 
model  with  n  =  2.1  and  po  =  3  x  10-12  g  cm-3  for  (3  Psc. 

4.3.  v  Cygnus 

v  Cyg  has  long  been  recognized  as  a  Be  star  with  strong  Ha 
emission  (Fleming  1891;  Campbell  1895).  Although  several  au¬ 
thors  have  noted  that  this  emission  stays  relatively  constant  over 
periods  as  long  as  decades  (see,  for  example,  Peters  1979),  others 
have  noted  periodic  outbursts  in  Ha  emission,  typical  of  many 
Be  stars  (Neiner  et  al.  2005).  The  interferometric  observations 


Spatial  Frequency  (tO6  cycles/radian) 

Fig.  7. — Interferometric  data  and  model  with  n  =  2.1  and  po  —  3.0  x  10  12  g 
cm-3  for  3  Psc.  The  dashed  line  represents  the  central  star,  and  the  solid  line  is  ob¬ 
tained  by  taking  a  Fourier  transform  of  a  synthetic  image  of  the  preferred  model. 
[See  the  electronic  edition  of  the  Journal  for  a  color  version  of  this  figure.] 
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Fig.  8. — Observed  Ha  line  for  v  Cyg  compared  with  the  best  predicted  lines  from  the  models,  which  are  also  in  agreement  with  the  interferometric  observations.  Solid 
lines  represent  the  model  profiles,  and  the  circles  represent  the  observed  line.  See  Table  3  for  the  model  parameters,  Ha  equivalent  widths,  and  the  reduced  \2  values  from 
interferometry  corresponding  to  these  models.  [See  the  electronic  edition  of  the  Journal  for  a  color  version  of  this  figure .] 


for  v  Cygnus  are  from  12  nights  between  2005  August  28  and 
September  26,  and  the  Ha  line  profile  was  acquired  on  2005 
September  1 6  (see  Table  1 ).  We  note  that  we  also  have  observa¬ 
tions  from  2005  October  1 1  that  confirm  the  spectroscopic  stabil¬ 
ity  of  the  star. 

Neiner  et  al.  (2005)  suggested  that  the  outbursts,  which  occur 
over  periods  of  a  few  years,  are  the  result  of  multiperiodic,  non- 
radial  pulsations.  Neiner  et  al.  (2005)  derived  the  stellar  param¬ 
eters  for  this  star  by  detailed  modeling  constrained  by  photometric 
and  spectroscopic  data.  Three  of  their  models  include  veiling  ef¬ 
fects  and  gravitational  darkening  for  a  range  of  stellar  rotational 
velocities:  80%,  90%,  and  95%  of  the  critical  velocity  (see  Neiner 
et  al.  2005  for  more  detail).  We  have  chosen  to  adopt  the  stellar 
parameters  from  their  model  D,  which  has  an  intermediate  value 
of  the  rotational  velocity  (90%  of  critical  velocity).  Also,  Neiner 
et  al.  (2005)  find  a  value  for  the  inclination  angle  of  i  ~  27°.  The 
stellar  radius  listed  in  Table  2  corresponds  to  the  equatorial  ra¬ 
dius  given  by  Neiner  et  al.  (2005).  This  is  a  reasonable  choice, 
since  the  star  is  nearly  pole-on.  The  stellar  parameters  adopted  for 
v  Cyg  are  reasonable  compared  to  other  values  presented  in  the 
literature  of  this  spectral  type  (see,  for  example,  Cox  2000). 

The  procedure  for  selecting  the  grid  of  models  to  be  com¬ 
pared  with  interferometric  observations  was  obtained  by  the 
same  selection  process  used  for  the  two  other  stars  in  this  in¬ 
vestigation.  The  final  selection  of  parameters  to  be  searched 
ranged  from  n  =  1.8  to  4.5  and  from  po  =  1.0  x  10~12  to  4.0  x 
10~9  g  cm~3,  yielding  a  total  of  193  models  constructed  for 
comparison  with  interferometric  observations.  The  correspond¬ 


ing  x2  values  had  a  minimum  of  1 . 14,  corresponding  to  the  model 
with  n  =  4.0  and  po  =  8.0  x  10~10  g  cm~3,  and  a  maximum  of 
471,  corresponding  the  model  with  n  =  1.8  and  po  =  8.0  x 
HH2  g  cm~3.  The  x2  values  revealed  several  models  that 
statistically  fit  the  observations  equally  well.  Again,  we  turn  to 
spectroscopy  to  remove  the  degeneracy. 

In  order  to  predict  Ha  profiles,  we  adopt  i  =  30°,  which  is 
consistent  with  Neiner  et  al.  (2005).  We  calculate  the  Ha  equiv¬ 
alent  widths  for  a  subset  (^22  models)  chosen  based  on  the 
interferometric  reduced  x2  <  2.  Statistically,  the  models  with 
larger  x2  am  poorer  fits  to  the  interferometric  data.  The  equiv¬ 
alent  widths  for  the  subset  ranged  from  a  minimum  of  —19.0  A 
to  a  maximum  of  —43. 1  A.  The  observed  Ha  line  has  an  equiva¬ 
lent  width  of  —27.3  A,  and  there  were  three  models  with  reasonable 
line  shapes  and  equivalent  widths  that  were  close  to  the  observed 
value  and  had  statistically  similar  interferometric  x2  values. 
Figure  8  shows  the  Ha  profiles  corresponding  to  these  models. 

TABLE  3 

Best-Fit  Models  for  v  Cygnus 


po  Ha  Equivalent  Width 

n  (g  cm-3)  (A)  X2a 

4.0 .  8.0  x  10-10  —26.0  1.14 

4.1  .  1 .0  10  "  -24.4  1.16 

4.2  .  1.5  xlO-9  -23.7  1.30 


3  Based  on  comparison  with  interferometry. 
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Fig.  9. — Observed  Ha  lines  for  v  Cyg  compared  with  the  best  predicted  lines  for  n  —  3.6  and  n  —  4.4.  Solid  lines  represent  the  model  profiles,  and  the  circles  represent 
the  observed  line.  The  model  on  the  left  has  a  lower  n  and  higher  density,  resulting  in  an  Ha  profile  that  is  too  strong.  The  model  on  the  right  has  a  higher  n  and  lower 
density,  resulting  in  an  Ha  profile  that  is  too  weak.  [See  the  electronic  edition  of  the  Journal  for  a  color  version  of  this  figure .] 


Spatial  Frequency  (108  cycles/radian) 

Fig.  10. — Interferometric  data  and  model  with  n  —  4.0  and  po  —  8.0xl0_1°g 
cm-3  for  v  Cyg.  The  dashed  line  represents  the  central  star,  and  the  solid  line  is 
obtained  by  taking  a  Fourier  transform  of  a  synthetic  image  of  the  preferred  model. 
[See  the  electronic  edition  of  the  Journal  for  a  color  version  of  this  figure .] 


Table  3  lists  the  values  of  n,  po,  Ha  equivalent  width,  and  the 
reduced  \2  based  on  comparisons  from  interferometry  for  these 
three  models.  The  models  with  /;  =  4. 1  and  n  =  4.2  in  Figure  8 
have  predicted  profiles  that  are  too  narrow  and  also  do  not  fit  the 
peak  of  the  observed  profile.  Since  the  \2  values  are  statistically 
similar,  we  adopt  the  model  corresponding  to  n  =  4.0  and  po  — 
8.0  x  10-1°  g  cm-3  in  Figure  8  as  the  best  fit.  As  a  further  illus¬ 
tration,  Figure  9  shows  Ha  profiles  corresponding  to  models  with 
smaller  and  larger  n  than  those  given  in  Table  3.  Specifically, 
these  profiles  correspond  to  models  with  («  =  3.6,  po  =  3.0  x 
10-1°  g  cnr3)  and  ( n  =  4.4,  po  =  2.0  x  10-9  g  cm-3),  with  re¬ 
duced  x2  values  of  1.37  and  1.29,  respectively.  The  shapes  of 
these  profiles  do  not  match  the  observed  line.  The  predicted 
profile  corresponding  to  n  =  3.6  has  too  much  emission  (with  an 
equivalent  width  of  —33.7  A)  and  is  too  wide  in  the  line  wings. 
The  predicted  profile  corresponding  to  n  =  4.4  has  too  little  emis¬ 
sion  (with  an  equivalent  width  of  —20.3  A),  and  much  of  the 
profile  is  narrower  than  the  observation.  Figure  10  shows  a  com¬ 
parison  of  the  interferometric  observations  and  the  best-fit  model 
for  v  Cyg. 

5.  CONCLUSIONS 

In  this  work,  disk  density  models  were  obtained  for  the  clas¬ 
sical  Be  stars  k  Dra,  /3  Psc,  and  v  Cyg  by  matching  the  observed 
interferometric  Ha  visibilities  with  Fourier  transforms  of  syn¬ 
thetic  images  produced  by  theoretical  models.  It  was  demon¬ 
strated  that  the  additional  constraint  obtained  by  the  (cotemporal) 
observed  Ha  line  profile  was  critical  to  selecting  a  density 
model.  This  was  particularly  evident  in  the  case  of  k  Dra.  We 
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TABLE  4 

Summary  of  Best-Fit  Models 


Po 

Star 

n 

(g  cm-3) 

k,  Dra . 

4.2 

1.5  x  10-10 

/3  Psc . 

2.1 

3.0  xlO-'2 

v  Cyg . 

4.0 

8.0  xlO-10 

have  demonstrated  that  in  the  case  of  these  three  stars,  the  de¬ 
generacy  in  the  interferometric  best-fitting  models  can  only  be 
removed  by  additional  observational  constraints. 

The  best-fit  models  to  both  the  interferometric  observations 
and  the  observed  Ha  line  profiles  for  the  three  stars  are  sum¬ 
marized  in  Table  4.  The  range  of  n  is  consistent  with  previous 
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